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FOREWORD 

This  collection  consists  of  papers  initially 
presented  by  the  authors  at  IEEE  Oceans  Conferences 
held  in  1972,  '73,  '74,  and  '76.  It  traces  the  de- 
velopment of  their  work  on  the  low  frequency  electro- 
magnetic signatures  of  perfectly  conducting  objects 
of  arbitrary  but  smooth  shape  immersed  in  a lossy 
medium;  i.e.,  the  ocean.  The  papers  summarize  their 
successful  effort  in  applying  physical  optics  approxi- 
mation techniques  to  a broad  class  of  problems  and  in 
validating  this  approach,  as  well  as  the  limits  of 
its  applicability,  by  comparison  with  the  exact  solu- 
tion for  a spherical  scatterer.  This  work  forms  the 
background  for  current  efforts  to  extend  this  approach 
to  scatterers  with  edges,  where  the  radii  of  curvature 
are  small  compared  to  the  skin  depth  in  the  medium. 
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ELECTRIC  AND  MAGNETIC  DIPOLE 
SCATTERING  FROM  A PERFECTLY  CONDUCTING 
SPHERE  IN  A HOMOGENEOUS 
LOSSY  MEDIUM 

D.M.  Bolle 
R.A.  Bowden 
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Abstract 

A comparison  of  the  use  of  magnetic  or 
electric  dipole  sources  in  geophysical  ex- 
ploration and  the  detection  of  conducting 
objects  in  a lossy  environment  such  as  the 
ocean  was  made  through  numerical  evaluation 
of  signal  levels  for  orthogonally  position- 
ed electric  and  magnetic  transmitter-re- 
ceiver systems  traveling  along  linear 
paths  near  such  objects. 


Introduction : The  problem  of  detection  of 

spherical  bodies  in  a lossy  environment  has 
received  considerable  attention.  An  early 
analysis  was  presented  by  March (D  in  1953. 
Wait  (2)  presented  a complete  analysis  of 
the  magnetic  dipole  in  the  proximity  of  a 
sphere  of  conductivity  and  permeability 
different  from  the  conducting  medium  in 
which  it  is  immersed.  Extensive  review 
articles  were  published  by  Ward  in  1967 
where  the  geophysical  application  was 
stressed'3,4'.  However  a dearth  of  pub- 
lished numerical  data  exists  except  for 
some  presented  in  ref.  6 on  plane  wave 
scattering  from  a sphere  in  a layered 
earth  and,  in  addition,  the  orthogonal  mag- 
netic dipole  transmitter/receiver  response 
on  the  earth's  surface  over  a buried  spher- 
ical conducting  sphere. 

The  work  reported  on  here  was  performed  to 
evaluate  the  relativemerits  of  electric  and 
magnetic  dipole  transmitter/receiver  sys- 
tems for  the  detection  of  spherical  con- 
ducting bodies  in  a lossy  medium.  The  nu- 
merical results  were  obtained  for  an  homo- 
geneous ocean  environment.  All  calcula- 
tions were  based  on  the  analysis  as  pre- 
sented by  D.S.  Jones'3'  which  follows  the 
earlier  work  of  C.T.  Tai  which  has  been  re- 
published recently. 

For  the  sake  of  completeness  the  basic  re- 
sults are  stated  for  the  magnetic  and  elec- 
tric dipole  sources. 


Analytical  Results:  The  electric  and  mag- 

netic  fields  for  the  dipoles  can  be  ex- 
pressed in  terms  of  the  electric  and  mag- 
netic Herzian  vector  potentials,  i.e., 

£ 2 W-77  +‘A*’/r« 

H x vx  TTC 

and  for  the  magnetic  dipole: 

£ = - iw^  7x  77^ 

£•  + klJTn 


where 


tlut-irj 

A*  C 

-y/rtf  ■ r 

~m  jfr  ' r 


and  where  r is  the  radial  coordinate  of 
a system  centered  on  the  dipole  source  and 
Ap  and  Bm  are  the  electric  and  mag- 
netic dipolemoments , respectively.  The 
lossy  nature  of  the  medium  is  included  in 
the  dielectric  constant;  thus 

c-  e '(/- V ir/wtf'J 

where  a and  e ' are  the  conductivity 
and  dielectric  constant  of  the  medium.  The 
propagation  constant  k is  given  by 

K -i 

or  , 

k » (>-'<■)/: 

where  5 = is  the  skin-depth  at  the 

frequency  a)  uS!nd  u is  the  permeability 
of  the  medium;  also,  we  have  used  the  fact 
a 


that 


we 


>>  1 for  the  lossv  medium  at 


the  frequencies  of  interest. 

The  fields  can  be  expressed  directly  ir. 
terms  of  source  functions  through  the  use 
of  tensor  Green's  functions'3'  . Thus,  for 
the  electric  dipole  (see  figure  I) , 
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E'-k'i 

h*  -iuA  e,.vxrl(#,J£) 

and  for  the  magnetic  dipole, 

c » - «'«/t  B rna  • vx  FL  (X'j  %) 

H,  £/*:«! 

where 

c .jr£L±li!:£**g& 

where , in  turn  “ “ 

for  R'  > R:  r =4  and  s = 1; 


and  for 


R'  < R:  r = 1 and  s = 4 , 


('>-/*•/)  !/[*  («+>)  ( »+H)l] 

KMC* 

C (Mi  % •)  »*  4 iK'J- 

1 

KM .j.M  . KM-  CM 

the  prime  on  p and  p indicates  that  r, 

9 and  $ are  replaced  by  r',6',  and 
It  is  readily  shown  that  £2  takes  the 
3ame  form  as  £^. 

For  either  type  of  dipole,  the  Green's 
tensor  consists  of  the  sum  of  two  parts; 
one  part  for  the  incident  field,  and  one 
for  the  scattered  field.  Thus, 

II  * T-il  + 4l 

and 

V-2  " ri2  + — s2 

where 

-SI  • ifKk  »»/  m--* 

f . Hr)  It)  I V.  J 

and  a and  3 are  unknown  coeffi- 
mn  mn 

cients.  r has  the  same  form  but  dif- 

ferent  coefficients  a__  and  b . 

mn  mn 


The  boundary  conditions , that  the  tangen- 
tial electric  field  must  vanish  on  the 
sphere,  are  given  by 

a x E = 0 
R 


(SR  * le  = ° 

sR  x (7  , rm)  = 0 


R = a.  This  leads  to  the  result 


* - J 

°v»/>  - “ 


where  Z = ka. 
Using  the  relation 


the  series  for  £.  can  be  summed  over  m 
to  obtain 

+ g(**v)l*'*v)K  (&.*') 

where,  for  R < R' 


CM  n /•“’') 
= 5 S'.  Cm] 

♦fj  n (hti)u 


2 


and 


y = a / d 
' n mn'  nrn 

5 = 3 /d 

n ran  mn 

When  R > R’ , R and  R'  are  inter- 
changed in  the  expressions  for  and  Vj. 

The  only  changes  for  the  magnetic  dipole 
case  are  in  and  where 


a /d 
mn  mn 


D 


n 


b / d 
mn 


mn 


replace 


and 


This  completes  the  formal  statement  of  the 
solution  to  the  problem  of  an  electric  or 
magnetic  dipole  in  the  presence  of  a con- 
ducting sphere  immersed  in  a homogeneous 
lossy  medium. 

Results:  All  computations  were  performed 

for  a medium  of  conductivity  of  4 mhos/ 
meter,  i.e.,  approximating  that  of  sea 
water.  The  geometric  dimensions  are  all 
normalized  with  respect  to  skin  depth 
which  is  used  as  the  characterizing  di- 
mension, while  the  field  strengths  are 
normalized  with  respect  to  the  dipole 
strengths  which,  for  the  electric  dipole, 
becomes 

— 5 3 E and  — 52H 

p - p 

while  for  the  magnetic  dipole  case  we  have 

— o3H  and  — 5 4E 

m — m — 


with  the  resulting  normalized  field 
strengths  expressed  in  db’s. 

Seme  sample  results  for  the  system  shown 
in  figure  2 are  given  in  figures  3(a)-'g) 
for  the  electric  dipole  source  with  an 
orthogonally  mounted  magnetic  dipole  sen- 
sor for  various  combinations  of  source- 
sensor separation  b and  sphere  radii  a. 

The  observed  signal  level  in  the  sensor  is 
shown  along  linear  paths  at  various  stand- 
off distances  s from  the  sphere. 

Figures  4 (a) -(g)  cover  the  same  situations 
but  with  a dipole  source  substituted  for 
the  electric  dipole.  In  appropriate  cases, 
i.e.  when  the  sphere  radius  is  large  com- 
pared to  the  baseline  or  source-sensor 
separation,  the  comparison  with  the  infi- 
nite plane  result  is  shown  for  reference 
purposes . 

For  the  first  case  the  results  conform  to 
expectations.  As  the  sphere  is  approached 
the  signal  strength  increases  monotonical- 
ly,  maximizes  at  the  nearest  point,  then 


monotonically  decreases  as  the  sphere  re- 
cedes. Only  when  the  baselengths  becomes 
of  the  order  of  the  skin  depth  and  the 
stand-off  distance  s is  small  do  depart- 
ures occur  from  this  simple  behavior  and 
the  curves  show  marked  asymmetry.  This 
effect  is  to  be  ascribed  to  the  fact  that 
there  occurs  strong  interaction  first  be- 
tween receiver  and  sphere  upon  approach 
followed  by  strong  interaction  between 
transmitter  and  sphere  as  one  recedes  which 
effects  are  separated  by  the  long  baseline 
of  the  system. 

In  the  second  case  it  is  seen  that  the 
signal  level  again  increases  monotonically 
upon  approach  but  at  a position  close  to 
the  nearest  point  the  spinal  level  decreas- 
es abruptly  to  a local  minimim.  This 
effect  is  only  observed  for  the  case  where 
the  source-sensor  separation  b is  much 
smaller  than  the  sphere  radius.  As  the 
baselength  b becomes  comparable  to  the 
sphere  radius  this  effect  is  no  longer  ob- 
served while  for  long  baselengths,  i.e., 
b v $ , one  obtains  approximately  Gaussian 
response  curves  again  but  for  the  small 
departures  at  small  stand-offs  s due  to 
the  effect  mentioned  earlier.  This  dip  in 
the  response  curve  for  small  baselengths 
to  sphere  radius  ratios  is  most  useful  for 
accurate  location  of  the  near  approach 
point  to  a highly  conducting  inclusion  in 
a conducting  medium  and  has  been  observed 
earlier  in  some  results  presented  in  refer- 
ence 4.  The  physical  basis  for  this  ef- 
fect is  understood  readily  upon  consulting 
figure  5 and  utilizing  the  image  dipoie 
approximation.  At  the  position  of  nearest 
approach  (position  A)  the  image  generates 
a relatively  small  signal  in  the  receiver 
but  as  the  systems  moves  away  from  this 
point  (position  B)  the  image  rotates  due 
to  the  sphere  curvature  with  a. resultant 
increase  in  the  received  signal  level  in 
the  orthogonally  oriented  sensor.  Of 
course,  as  the  system  moves  further  out, 
the  field  and  thus  the  received  signal  must 
eventually  become  weaker. 

Conclusion : Numerically  computations  were 
made  comparing  the  effectiveness  of  elec- 
tric and  magnetic  dipole  sources  in  source- 
sensor systems  for  the  detection  of  highly 
conducting  spherical  inclusions  in  a mod- 
erately conducting  environment  such  as  the 
ocean.  Differences  in  response  ware  noted 
for  the  magnetic  versus  electric  dipole 
transmitter  systems  where  the  source  and 
sensor  are  oriented  orthogonally  so  as  to 
yield  small  signals  in  a homogeneous  envi- 
ornment.  The  particular  feature  noted  was 
the  nearest  point  local  minimum  obtained 
for  the  magnetic  dioole  source  system. 
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ABSTRACT 

The  utility  of  the  physical  optics  approx- 
imation is  explored  for  the  computation  of  elec- 
tromagnetic fields  of  dipole  sources  in  the 
neighborhood  of  conducting  bodies  of  diverse 
shape  in  a lossy  environment  such  as  the  ocean. 
Excellent  agreement  with  previous  exact  results 
is  obtained  for  spheres  down  to  one  skin  depth 
diameter  with  source-detector  separation  up  to 
one  skin  depth.  Only  when  the  sphere  size  is 
reduced  to  one- tenth  skin  depth  do  errors  in 
excess  of  three  decibels  begin  to  apoear.  Such 
comparisons  have  confirmed  the  accuracy  of  the 
method  for  application  to  non-spherical  body 
shapes.  The  computer  program  allows  investigation 
of  the  response  for  orthogonally  oriented  source 
and  detection  dipoles  in  the  presence  of  smooth, 
curved,  conducting  bodies  in  a homogeneous 
lossy  medium.  Results  are  given  for  spheres, 
infinite  cylinders,  prolate  spheroids  and  super- 
elliptic  spheroids  exposed  to  magnetic  dipole 
sources . 

INTRODUCTION 

The  detection  and  identification  of  sub- 
merged bodies,  now  considered  a classical  prob- 
lem, has  seen  a great  deal  of  ingenuitv  applied 
to  theoretical  and  experimental  methods  for  its 
solution.  Numerous  systems  have  been  oroposed, 
many  have  seen  realization  (at  least  in  pilot 
form)  and  have  been  subjected  to  considerable 
v"-ement  and  sophistication.  Sonar  technology 
.as  .-een  heavily  dominant  in  this  field  but  in 
recent  years  the  electromagnetic  spectrum  has 
been  closely  scrutinized  for  possible  contribu- 
tions to  this  Comdex  field.  The  frecuencv 
bands  considered  range  from  the  optical  region 
(laser  technology)  to  low  frepuencv  elect omae- 
netic  waves.  Effective  utilization  of  the  far 
greater  part  of  the  electromagnetic  spectrum  is 
blocked  by  the  characteristic  nature  of  the 
medium.  The  scattering  and  loss  orooerties  of 
the  ocean  severely  limit  the  usable  freouen^  y 
bands.  This  meant  that  much  effort  was  direct- 
ed towards  Electrostatic  and  Maenetostatic 
modes  of  ooeration.  Both  sources  and  detectors 
have  been,  and  continue  to  be,  the  subject  of 
intensive  investigations. 


There  are,  in  addition  to  the  advantages 
of  using  static  electromagnetic  systems,  pecul- 
iar advantages  to  using  low  freouencv  time 
varying  electromagnetic  systems . One  obvious 
attractive  feature  is  the  inherent  advantage 
of  detection,  discrimination  and  processing  of 
time  varying  signals  in  a noisy  environment 
where  the  freouencv  is  precisely  known  and 
controlled.  In  addition,  as  will  be  shown  la 
this  paper,  it  is  possible  to  develop  signal 
signatures  for  a large  class  of  electrically 
highly  conducting  bodies  where  such  bodies  are 
submerged  in  the  ocean.  Contrary  to  the 
extreme  mathematical  and  computational  complex- 
ities which  arise  when  the  acoustic  backscatter- 
ed  signal  from  a submerged  bodv  of  complex 
shaped  must  be  obtained  it  will  be  show  that  a 
relatively  simple  methodology  suffices  when 
such  targets  are  considered  in  the  low  fre- 
auencv  electromagnetic  field. 

The  methodology  that  Is  being  developed 
here  ha6  the  Inherent  capability  of  supplement- 
ing current  standard  technioues  in  the  detec- 
tion and  identification  of  bodies  submerged  in 
shallow  as  well  as  the  deep  ocean  environment. 
Geophysical  exoloration , off-shore  technology 
and  many  other  associated  fields  can  benefit 
from  a detailed  knowledge  of  the  electromagnetic 
signatures  of  bodies  whose  conductivities  differ 
substantially  from  that  of  the  ocean  and  the 
ocean  bottom.  Search  and  rescue  ooerations  and 
the  recovery  of  lost  eauioment  and  material  in 
ocean  ooerations  would  stand  to  benefit  from 
the  further  developments  of  this  aspect  of  the 
technology. 

An  earlier  paoer  (Ref.  1)  presented  an 
exact  solution  for  spherical  conducting  bodies 
in  the  ocean.  Even  for  such  a simole  case  the 
computational  effort  was  considerable  and  a 
-ore  general  and  manageable,  though  approximate, 
■ethod  is  sought  for  bodies  of  more  complex 
shape . 

Ehvsicai  optics  suggested  itself  for  two 
reasons.  First,  in  a lossv  medium  the  influ- 
ence of  the  source  is  strong  onlv  over  a limit- 
ed part  (about  one  skin  depth  in  radius)  of  the 
illurinated  region.  Rurther,  anv  creeping  waves 
are  sufficiently  attenuated  so  that,  for  the 
class  of  smooth  scatterers  considered  here. 
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induced  surface  currents  may  be  reasonably 
approximated  by  physical  optics.  Second,  the 
physical  optics  approximation  is  simple  and 
readily  implemented  for  objects  of  arbitrary 
shape.  The  induced  surface  current  is  given 
by 

J * 2n  x Hi 


(1  ♦ ikr'  - k‘r'2)sin8'e"*,<r 


(r'/«r 


where 

I m | * magnetic  dipole  moment  (absolute  value). 


where  n is  the  outward  normal  and  Hj  is  the 
incident  magnetic  field  intensity.”  That  is,  it 
is  assumed  that  the  surface  current  induced  at 
a point  on  a curved  surface  is  that  which  would 
be  induced  on  an  infinite  flat  plane  tangent 
to  the  surface  at  that  point.  Thus,  physical 
optics  results  become  exact  for  any  geometry 
for  which  image  theory  will  also  yield  a 
solution. 


and  the  primed  variables  refer  to  the  spherical 
coordinate  system  centered  on  the  dipole  loca- 
tion with  the  z - axis  in  the  direction  of  the 
dipole.  The  propagation  factor  k is  given  by 

k = (1  - D/S 

1/2 

where  S * (2/wuo)  is  the  skin  depth  and 
ui  * angular  freouency 


Computer  subroutines  covered  geometries- 
for  the  infinite  plane,  infinite  cylinder, 
sphere,  oblate  and  prolate  spheroids,  and  a 
superelliptic  spheroid  with  an  elliptic  section 
on  the  longitudinal  axis  and  a superelliptic  - 
section  in  the  transverse  plane.  Excluded  at 
present  are  scatterers  with  sharp  edges  or 
ridges. 

Since  the  scattered  field  intensities  are 
given  in  terms  of  an  integration  over  the  in- 
duced surface  current  distribution,  the  surface 
must  be  suitably  scanned  and  the  integration 
replaced  bv  a summation.  The  integration  over 
the  surface  was  performed  by  projecting  points 
onto  the  scattering  surface  from  suitably  placed 
points  on  a plane  tangent  to  the  scattering 
surface  at  the  point  nearest  the  source.  Two 
checks  are  made:  one  to  ensure  that  only  the 
illuminated  part  of  the  surface  is  included  in 
the  summation,  and  a second  to  terminate  the 
process  when  the  contribution  from  the  last 
scan  circle  is  small  compared  to  the  total  sum. 

The  program  was  first  checked  for  the  in- 
finite plane.  For  one-tenth  skin  depth  incre- 
ments in  scan  circle  radii  and  in  scan  point 
spacing  on  the  scan  circles , agreement  with 
image  theory  was  within  one  percent. 

A comparison  is  also  made  between  the  phys- 
ical optics  approximation  and  the  exact  solution 
for  the  sphere.  Agreement  exceeded  expectation, 
and  noticeable  quantitative  error  occurs  only 
for  sphere  radii  of  one-tenth  skin  depth  or 
smaller  in  combination  with  transmitter-receiv- 
er base  lengths  of  one-tenth  skin  depth  or  larg- 
er and  even  in  these  cases  the  qualitative 
results  were  adeouate. 


THEORY 


The  field  components  of  a magnetic  dipole 
are  given,  for  a lossy  medium,  by, 

'ml  e -{kr' 

r 3 J-  • = r (1  ♦ ikr'  )sin8'e  1Kr 

io(r'/4K 


H ,«  — • — T (1  ♦ ikr,)cose'e'ikr' 

r 4 nSJ  (r ’/«) 


u * permeability  of  the  medium 
o * conductivity  of  the  medium 

The  field  intensities  are  exnressed  in  this  way 
so  that  the  magnetic  fields  components  H , and 
H, , may  be  normalized  with  respect  to  the 
auantity  |m|/4it63  . 

Having  obtained  the  induced  current 
J(  x.  ,yjj  ,z,  ) at  some  point  on  the  surface  of  the 
scattering  body,  and  the  area  element  ds,  the 
scattered  field  H (x,,yf,z.)  at  a given  point 
is  given  by 

, f -ikr 

H Z -i.  I J x 7 (e  /r  )ds 
-»  u*  ls  ~ ° 

where  r is  the  distance  separating  these  two 
points.  Or,  in  terms  of  a finite  summation 
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where  a is  the  unit  vector  directed  from  the 
o 

body  point  to  the  field  point  and  p,q  are  scan 
oarameters . 

Some  results  are  illustrated  in  figures  1 
through  6.  In  every  case  the  calculations  were 
performed  for  a system  using  a magnetic  dipole 
as  the  source  and  a magnetic  field  sensor 
oriented  normal  to  it  so  that  a zero  output 
results  in  the  absence  of  the  scattering  body. 
The  abscissa  shows  the  position  of  the  midpoint 
of  the  transmitter-receiver  baseline  as  a 
straight  line  path  is  traversed  past  the  per- 
turbing object.  The  parameter  associated  with 
each  curve  gives  the  distance  of  closest  ap- 
proach or  standoff.  The  ordinate  gives  the 
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normalized  magnetic  field  intensity.  The  solu- 
tion for  an  infinite  plane  scatterer,  when  shown 
for  comparison,  is  indicated  by  the  diamond 
shape . 

Figures  1 through  4 compare  the  physical 
optics  approximation  (circle  points)  with  the 
closed  form  solution  (solid  curves).  Figure  1 
shows  the  excellent  agreement  for  a sphere 
radius  a * 1.0,  baseline  b = 0.1  and  standoffs 
S from  0.1  to  1.0.  Figure  2 gives  results  when 
the  baseline  is  increased  to  1.0;  agreement  is 
still  excellent  although  a small  deviation  is 
observed  about  the  point  of  nearest  approach. 

It  may  be  noted  here  that  the  physical 
optics  results  may  be  as  accurate  as  the 
closed  form  solution.  In  the  latter  case  con- 
vergence problems  arise  when  source  and  re- 
ceiver are  the  same  distance  from  the  center 
of  the  sphere ; i . e . , in  the  present  case , when 
the  system  is  at  or  near  the  point  of  closest 
approach . 

Figures  3 and  4 illustrate  the  effect  of 
reducing  the  sphere  radius  to  0.1.  It  may  be 
observed  that  agreement  with  the  exact  results 
is  reasonable  around  the  nearest  approach  for 
small  standoff  and  the  error  increases  both 
for  greater  standoff  and  when  the  system  is 
furthest  removed  from  the  point  of  nearest 
approach.  The  maximum  error  is  at  most  10  dB 
for  a signal  strength  of  50  dB  below  peak 
response. 

It  may  also  be  noted  in  connection  with 
Figure  3 that  several  scanning  mesh  parameters 
were  tried  varying  from  .1  to  .01  skin  depth. 

A noticeable  change  was  observed  in  going  from 
0.1  down  to  0.03  but  beyond  that  any  further 
mesh  refinement  has  negligible  effect. 

The  curves  of  Figure  5 give  physical 
optics  results  for  the  infinite  cylinder  of  unit 
radius  with  the  system  traverse  occurring  per- 
pendicular to  its  axis.  Results  are  very  simi- 
lar to  that  for  the  sphere.  This  illustrates 
the  limited  angular  region  over  which  the  mag- 
netic dipole  is  effective  and  confirms  earlier 
work  reported  in  Reference  2. 

Figure  6 gives  results  for  the  super- 
ellipsoid:  ellipsoidal  in  longitudinal  section 
and  super-elliptic  of  order  4 in  the  transverse 
section.  The  resulting  curves  have  the  form  to 
be  expected  for  a scatterer  which  is  broader 
and  flatter  than  the  sphere.  It  should  also  be 
realized  that  the  "corner"  of  the  super-ellipse 
may  not  be  adequately  represented  by  a scan 
parameter  of  one-tenth  skin  depth,  resulting 
in  some  smoothing  of  the  response  curve. 

CONCLUSIONS 

(i)  The  physical  optics  approximation  gives 
excellent  results  over  a wide  range  of  para- 
meters when  applied  to  scattering  in  a lossy 
medium. 

(ii)  The  range  of  parameters  over  which 
this  theory  may  be  used  coincides  with  those 


of  interest  in  the  studv  of  conducting  inhomo- 
geneities in  a lossy  environment. 

(iii)  Both  qualitatively  and  quantitatively , 
the  results  for  the  sphere  may  be  extrapolated 
to  a large  degree  to  yield  information  on  more 
complex  shapes  such  as  the  spheroid  with  super- 
elliptic  transverse  cross  section. 

(iv)  A method  has  been  investigated  and 
proven  accurate  for  predicting  the  electromag- 
netic signature  of  conducting  objects  in  the 
ocean  environment  for  object  sizes  down  to  a 
fraction  of  the  skin  depth  for  the  given 
signal  freouencv. 

(v)  The  nature  of  the  signatures  obtained 
shows  that  accurate  position  indication  may 
be  obtained  using  this  system. 
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INTRODUCTION 

In  previous  papers  (1,2)  th«  utility  of  the 
physical  optics  approximation  in  predicting  the 
electromagnetic  signatures  of  conducting  bodies  im- 
mersed in  the  ocean  was  investigated.  It  has  been 
shown  that  the  physical  optics  approximation  gave 
excellent  agreement  with  exact  results  for  spheres 
or  radius  of  one  skindepth  or  larger  of  the  elec- 
tromagnetic radiation  in  sea  water  and  good  results 
for  spheres  down  to  a radius  of  one  tenth  of  a skin- 
depth.  The  greatest  disagreement  incurred  between 
the  exact  and  approximate  results  was  about  three 
dS  in  the  latter  case. 

The  results  were  subsequently  extended  to  more 
complex  bodies  such  as  the  cylinder,  ellipsoid  and 
superellipsoid. 

In  this  work  the  results  obtained  in  the  next 
logical  sequential  step  in  generalizing  this  ap- 
proach to  objects  with  sharp  edges  will  be 
presented. 

DISCUSSION 

It  is  desirable  to  be  able  to  predict  the  sig- 
nature, not  only  of  smooth  objects,  but  also  of 
composite  objects  which  are  largely  smooth  but  may 
have  sharp  edges  and  sharp,  acute-angled 
protrusions. 

To  examine  the  utility  of  the  physical  optics 
approximation  when  the  radius  of  curvature  of  parts 
of  the  surface  of  the  objects  are  much  smaller  than 
the  skindepth,  the  magnetic  dipole  in  the  vicinity 
of  the  edge  of  a semi- infinite  conducting  plane 
was  considered. 

This  geometry  can  serve  as  a canonical  problem 
and  produces  the  most  severe  test  for  the  physical 
optics  approximation.  Secondly,  exact  solutions 
are  available  (3,u)  and  thus  the  necessary  compar- 
ison can  be  made.  It  is  expected  that  the  physical 
optics  approximation  will  not  suffice  for  all  ori- 
entations of  the  dipole  with  respect  to  the  direc- 
tion of  the  edge  of  the  half-plane.  Thus,  the 
availability  of  exact  solutions  will  allow  the 
evalutations  of  a modified  physical  optics  approach. 


the  master  program  and  calculate  directly  the  re- 
sponse of  a large  class  of  objects  exposed  to  a 
magnetic  dipole  source  and  an  orthogonally  oriented 
magnetic  field  component  sensor. 

THEORY 

The  physical  optics  approximation  consists  of 
replacing  the  exact  current  distribution  at  a point 
on  a curved  conducted  body  by  the  surface  current 
that  would  be  induced  on  a infinite  plane  if  it  were 
tangent  to  the  curved  body  at  the  point  of  interest. 
Thus,  for  an  infinite  plane  the  result  is  exact  and, 
as  the  curvature  increases,  we  can  expect  deterior- 
ation in  the  results  obtained.  The  fact  that  we 
are  considering  a system  immersed  in  a lossy  medium, 
i.e.,  sea  water,  accounts  for  the  fact  that  we  ob- 
tain excellent  to  good  results  down  to  large  curva- 
tures, i.e.,  radii  of  curvature  down  ti  one  tenth 
of  a skindepth.  At  a sharp  edge  we  must  expect 
gross  errors  since  it  is  well-known  that  surface 
current  component  parallel  to  an  edge  should  ex- 
hibit singular  behavior  and,  for  a semi-finite 
plane , the  current  component  parallel  to  the  edge 
will  increase  as  r*^-^2  where  r is  the  distance  from 
the  edge.  Surface  current  components  normal  to  the 
edge  must  vanish  as  r goes  to  zero.  However, 
physical  optics  predicts  constant  values  for  both 
components  in  the  vicinity  of  an  edge. 

To  construct  a Modified  Physical  Optics  approx- 
imation presents  no  great  difficulty.  It  is  an 
easy  matter  to  write  suitable  approximations  for 
the  current  components  in  terms  of  two  parameters 
then  four,  etc.,  until  we  have  achieved  an  iccuracy 
adequate  to  our  purpose. 

Previous  studies  of  surface  currents  near  the 
edge  of  acute  angled  wedge-shaped  semi-infinite  con- 
ductors illuminated  by  a plane  electromametic  wave 
give  considerable  insight  into  the  approximations 
to  be  employed  ( 5 ) . 

Thus,  the  physical  optics  approximation 

J = 2n  x (1) 

where  J is  the  surface  current  density,  n the  sur- 
face normal  and  H.  the  incident  magnetic  field 
intensity,  is  modified  to  become 


Once  an  adequate  model  has  been  established  j a 2n(n  x H.)  • (t  x n)(l  - e"or)(n  x t) 

for  predicting  the  current  distributions  induced  — -— i - - --  (.) 

near  edoes.  we  can  introduce  anrh  model  lino  into  . - / _ u \ . _/ , . ,_-l/2,_ 
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where  r Is  a unit  vector  parallel  to  the  edge.  As 
the  half-plane  is  allowed  to  assume  a wedge  shape 
so  a and  A will  change  as  will  the  exponent  of  r in 
the  second  term  for  J. 

RESULTS 

Computer  programs  were  prepared  yielding  the 
magnetic  field  sensed  by  a detector  oriented  normal 
to  the  magnetic  dipole  at  a fixed  distance  from  the 
dipole  source  as  this  source  traversed  a path  normal 
to  the  edge  of  the  semi-infinite  conducting  plane. 

The  methods  employed  were;  first,  the  physical  optics 
approximation;  secondly,  exact  results  based  on 
work  by  3owman  and  Senior  (6);  and,  thirdly,  the 
modified  physical  optics  approximation  as  given 
by  equation  (2). 

The  two  figures  give  samples  of  the  computed 
results.  The  path  of  the  dipole  source  — sensor 
system  was  normal  to  the  edge  of  the  half-plane  and 
parallel  to  that  plane.  The  separation  between  the 
source  and  the  sensor  was  fixed  at  one  tenth  of  a 
skindepth  for  these  calculations  and  the  dipole 
oriented  normal  to  the  semi- infinite  plane 
(x  > 0,  y = 0). 

Figure  1 shows  sample  calculations  obtained 
from  the  exact  theory  for  that  component  of  the 
suitable  normalized  magnetic  field  intensity  which 
is  zero  in  the  absence  of  any  perturbing  object  in 
the  vicinity  of  the  dipole.  Above  the  semi- infinite 
plane  and  away  from  the  edge  the  results  rapidly 
approaches  the  limit  value  obtained  from  image 
theory  for  the  magnetic  dipole  above  an  infinite 
perfectly  conducting  plane,  while  beyond  the  edge 
the  results  again  rapidly  approach  the  asymptotic 
slope  of  30  dB's/skindepth.  The  'edge-effect'  is 
evidently  limited  to  a region  not  extending  beyond 
one  and  a half  times  the  stand-off,  i.e.,  the  dis- 
tance of  the  path  above  the  surface  y - 0,  which 
contains  the  semi- infinite  plane. 

Figure  2 shows  the  results  of  the  calculations 
based  on  the  exact  theory  (curve  A)  for  the  case 
where  the  stand-off  is  half  a skindepth  and  compares 
it  to  the  results  of  the  physical  optics  approxima- 
tion (curve  B)  as  well  as  the  results  obtained  from 
the  Modified  Physical  Optics  approximation  with  the 
parameters  values  A = 2.0,  a = .5,  (curve  C),  and 
A = .1,  a * 1.0  (curve  D).  It  is  seen  that  the  two 
parameter  approximations  allow  wide  latitude  in 
increasing  the  accuracy  of  the  physical  optics  ap- 
proximation and  models  the  behavior  of  the  system 
in  the  neighborhood  of  the  edge  appropriately. 

It  is  concluded  that  the  physical  optics  method 
may  be  modified  to  include  objects  with  sharp  edges 
and  that  this  method  is  therefore  sufficiently  flex- 
ible to  allow  computation  of  the  electromagnetic 
signature  of  complexly  shaped  perfectly  conducting 
metallic  objects  in  the  ocean. 
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Abstract 

The  objective  of  this  study  vas  to  develop 
an  approximate  method  for  predicting,  with 
acceptable  accuracy,  the  electromagnetic  field 
scattered  by  perfectly  conducting  bodies  of 
complex  shape  Immersed  in  a lossy  medium,  i.e. 
the  ocean,  when  subjected  to  the  field  of  an 
oscillating  magnetic  dipole  source. 


Introduction 

Since  the  class  of  scatterers  of  Interest 
had  characteristic  dimensions  of  the  order  of, 
or  larger' than,  the  electromagnetic  wavelength 
or  the  'akin  depth'  in  the  lossy  medium  and, 
furthermore,  were  initially  restricted  to  those 
having  smooth  bounding  surfaces,  it  was  decided 
to  investigate  the  physical  optics  approxima- 
tion. The  approach  appeared  particularly 
attractive  because  of  the  ready  adaption  of  the 
computer  program  to  handle  a wide  variety  of 
scattering  shapes. 

The  study  consisted  of  four  Interlocking 
parts.  First,  it  was  necessary  to  determine 
the  effectiveness  and  accuracy  of  the.  physical 
optics  approximation  in  this  environment.  To 
this  end  not  only  was  a computer  program 
developed  which  allowed  the  calculation  of  the 
scattered  field  from  a smooth,  but  as  yet  un- 
specled,  body  shape  but  also,  using  a formulation 
by  3.  S.  Jones  Cl],  the  exact  solution  for  the 
magnetic  dipole  field  scattered  by  a perfectly 
conducting  sphere  lsnersed  in  a conducting 
medium  was  programmed  and  obtained. 

Some  sample  results  obtained  with  this  pro- 
gram are  shown  in  figure  (1) . These  results  are 
taken  from  earlier  work  [2]  which  we  reported  on 
at  Ocean  '72.  The  results  shown  are  for  the 
case  where  we  have  a magnetic  dipole  source  and 
a magnetic  field  sensor  pair  passing  the  spher- 
ical inclusion  with  the  magnetic  dipole  oriented 
in  the  direction  of  motion  and  the  magnetic 
sensor  directed  normal  to  the  path  and  in  the 
plane  containing  the  coordinate  center  of  the 
sphere.  The  seperatlon  between  the  dipole  and 
sensor  is  .Oil  for  (la)  and  (lb)  and  .If  for 
(lc)  and  (Id)  while  the  minimum  approach 


distance  S is  shown  as  a parameter  in  all 
four  cases.  The  pertinenc  sphere  radii  are 
shown  in  the  figures,  and  5 is  the  skin 
depth  in  the  medium. 

The  second  phase  of  the  study  vas  to  ob- 
tain data  which  would  allow  us  to  evaluate 
the  range  of  the  parameters  for  which  the 
physical  optics  approximation  gave  results 
which  were  adequate,  say  within  three  db, 
when  compared  with  the  exact  results  obtained 
for  the  spherical  inclusions.  The  physical 
optics  approximation  assumes  chat  on  curved 
perfectly  conducting  bodies  the  total  local 
tangential  magnetic  field  intensity  is  cvlce 
that  of  the  incident  field.  Thus,  for  an 
Infinite  plane  Che  physical  optics  result  is 
exact  but  the  approximation  deteriorates  as 
the  curvature  of  the  surface  increases.  In 
this  sense  the  comparison  with  the  sphere  gives 
us  an  ideal  canonical  situation  in  which  we  can 
evaluate  the  effect  of  curvature  on  accuracy. 
Figure  (2)  gives  some  sample  results  which 
compare  the  physical  optics  method  with  Che 
exact  results. 

It  was  surprising  to  find  that  in  the  case 
when  the  conductivity  o of  the  medium  was 
taken  to  be  A mhos /meter  the  agreement  between 
the  approximate  and  exact  results  did  not  begin 
to  break  down  until  we  had  reached  a sphere 
size  of  only  .1  5 (see  reference  [3]). 

Armed  with  the  confidence  gained  from  the 
results,  data  were  obtained  for  a variety  of 
more  complex  surfaces  such  as  the  infinite 
cylinder  and  the  superelliptlc  spheroid  (see 
figures  (3)  and  (4)).  The  computer  program  was 
cast  in  such  a form  that  additional  geometries 
could  be  Introduced  readily  as  subroutines. 
Currently  available  are  (1)  the  infinite  plane, 
(2)  the  infinite  cylinder,  (3)  the  sphere,  (4) 
the  prolate  or  oblate  spheroid,  (5)  the  super- 
elliptic  spheroid,  (6)  Che  hyperbolic  cylinder. 

The  reason  for  Including  the  infinite  plane 
was  to  allow  the  numerical  procedures  to  be 
checked  under  the  condition  where  the  physical 
optics  approach  does,  in  fact,  yield  the  exact 
and  correct  result  thereby  isolating  procedural 
and  numerical  sources  of  error  from  those 
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introduced  due  to  surface  curvature  where  the 
P.  0.  approach  13  not  exact. 

The  third  part  of  the  research  program 
consisted  of  extending  the  physical  opcics 
approach  to  bodies  with  edges. 

A direct  extension,  once  a suitable  canon- 
ical geometry  for  testing  purposes  was  chosen, 
was  rather  straightforward  in  that  only  a 
further  subroutine  needed  to  be  introduced  into 
the  master  program.  The  choice  of  the  canonical 
geometry  was  limited  to  those  few  for  which 
exact  solutions  already  existed.  For  this  reason 
the  semi -infinite  perfectly-conducting  plane  was 
chosen.  Two  approaches  leading  to  exact  solu- 
tions are  described  in  the  literature  f4J. 

Some  initial  results  were  generated  which 
made  use  of  the  direct  physical  optics  method. 
Sample  results  are  shown  in  figures  (5)  and  (6) . 
It  is  clear,  however,  that  the  physical  optics 
approximation  does  not  generate  the  appropriate 
singularity  behavior  in  the  field  structure 
near  the  edge.  For  example,  in  the  case  of  the 
semi-infinite  plane  the  component  of  the  field 
normal  to  the  edge,  l.e.,  the  surface  current 
density  component  parallel  to  the  edge,  must 
behave  as  r"l/2  as  the  edge  is  approached  [5] 
whereas  the  toul  surface  current  density  nor- 
mal to  the  edge  must  vanish.  The  program  was 
therefore  modified  to  allow  the  introduction  of 
two  parameters  which  could  subsequently  be  ad- 
justed to  minimize  the  discrepancies  between 
the  exact  and  approximate  results  and  allow  the 
extension  of  the  P.  0.  approach  to  bodies  with 
corners  and  edges,  l.e.,  the  induced  surface 
current  Jg  was  given  the  form 

Jg  * 2(n  x Hj)  • {(t  x n)  (1  - e'ar)(n  x t) 

+ t(1  + Ar  ^*)t} 

where  n and  T are  unit  vectors  normal  and 
tangential  Co  the  edge  of  Che  half-plane,  B^ 
in  the  incident  magnetic  field  intensity,  a 
and  A are  che  two  adjustable  parameters. 

Figure  (7)  shows  the  effect  of  changes  in  these 
parameters. 

The  fourth  and  final  part -of  the  current 
program  consists  of  obtaining  numerical  results 
from  the  exact  solutions  for  the  canonical  geo- 
metry. This  effort  is  complicated  by  the  fact 
that,  as  in  che  case  for  che  exact  solution  for 
the  spherical  geometry,  special  functions  must 
be  generated  (Hankel  functions)  which  have  com- 
plex arguments.  At  this  writing  these  results 
have  Just  become  available.  In  figure  (8), 
however,  we  display  the  exact  solution  for  che 
semi-infinite  plane  but  for  the  edge  scattered 
term  which  has  not,  as  yet,  been  included. 

Conclusion 

This  paper  discusses  the  current  status  of 
an  investigation  into  che  use  of  the  physical 
optics  approximation  in  magnetic  dipole 
scattering  from  perfectly  conducting  inclusions 


in  an  imperfectly  conducting  medium.  The  in- 
clusions are  generally  taken  to  be  smooch  but 
may  In  addition  exhibit  edges  and  comers. 
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Figure  2:  The  Physical  Optics  Approximation  and  the  Exact  Solution  Compared. 

Spere  radius  - a.  S ■ nearest  approach  distance. 
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